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We consider the coupling of a qubit in a pure state to an environment in an arbitrary state, 
and characterize the possibility of qubit-environment entanglement generation during the evolution 
of the joint system, that leads to pure dephasing of the qubit. We give a simple necessary and 
sufficient condition on the initial density matrix of the environment together with the properties 
of the interaction, for appearance of qubit-environment entanglement. Any entanglement created 
turns out to be detectable by the Peres-Horodecki criterion. Furthermore, we show that for a large 
family of initial environmental states, the appearance of nonzero entanglement with the environment 
is necessarily accompanied by a change in the state of the environment (i.e. by the back-action of 
the qubit). 


When a quantum system, a qubit (Q) in the context of 
this paper, is coupled to an environment (E), an initial 
pure state of the qubit evolves into a mixed state. It is 
widely recognized that there is an intimate relation be¬ 
tween this process of decoherence Qi and creation of 
qubit-environment entanglement (QEE). A precise state¬ 
ment can be made for E initialized in a pure state: the 
establishment of QEE is then equivalent to the reduction 
of purity of the reduced density matrix of Q. However, 
in the case of an initial mixed state of E, the situation is 
more complicated: the state of Q can lose its purity while 
no QEE is established Q. This should not be surprising, 
since QEE should be associated with decoherence, un¬ 
derstood as a process in which information is transferred 
from Q to E (i.e. E is in some sense “measuring” Q), not 
with the bare fact that the state of Q is becoming mixed. 
The latter can happen when there is no influence of Q 
on E. For example, when the self-Hamiltonian H-e of E 
commutes with the qubit-environment interaction Vqe, 
and when the initial state of E fulfills pe(0) = /{He) 
(e.g. it is a thermal state determined by He), there is no 
back-action of Q on E, while E is simply a source of ran¬ 
dom, classical, and quasi-static fields acting on Q 
In this case of so-called random unitary (RU) evolution, 
the purity of Q decays while no QEE is established. 


A natural question to ask is whether the RU case is the 
only one for which QEE does not accompany the loss of 
purity of Q. While specific kinds of environments and 
qubit-environment couplings were investigated in this 
context (e.g. quantum Brownian motion [3, 0 or pure 
dephasing due to a bath of noninteracting bosons i), 
the general answer to this question seems to be lacking. 
This is to a large degree caused by the fact that quan¬ 
tification (or even checking for presence) of QEE is a 
very hard problem when mixed states of the total system 
are considered, and when the dimension of Hilbert space 
of E is larger than 3 Ifil-T^. The existence of bound 


entanglemet 14, 1^ (which is not detected by the Peres- 
Horodecki criterion 16|, ll7| of negativity of a partially 
transposed matrix of the total system) severely hampers 
the task of general understanding of QEE [l^ . 


In this paper we present a complete characterization of 
conditions for the appearance of QEE in a less general, 
but physically well-motivated situation. Specifically, we 
focus on the case of pure dephasing (PD) of the qubit, 
defined by condition that [Hq , Vqe] = 0 (where Hq is 
the self-Hamiltonian of Q). The PD case is not only the 
paradigmatic example for relation between decoherence 
of Q and establishment of QEE [l| , it is also often encoun¬ 
tered in experiments, when the energy splitting of Q is 
tuned to values for which the exchange of energy between 
Q and E is suppressed either because of diminished cou¬ 
pling, or because of diminished density of states of E with 
matching energies. In such a constrained setup we derive 
the necessary and sufficient condition for nonzero QEE. 
This condition is that for the initial state of the environ¬ 
ment /Oe(0) = Cn \n) (n|, QEE is absent at time t iff for 
any states |i) and |n) for which Cif^Cn the qubit-induced 
evolution evolution of E does not couple these states (see 
below for a more precise formulation). This condition 
shows the lack of QEE for a completely mixed pE (0) (this 
statement is not trivial, as one may think, since there ex¬ 
ist non-pure-dephasing evolutions that lead to QEE when 
the initial density matrix of the environment proportional 
to unity, see Appendix for an example), for the RU case, 
and also for a family of cases which correspond to a mix¬ 
ture of the former two (i.e. when both pe(0) and Vqe 
are nontrivially constrained). Furthermore, QEE is al¬ 
ways detected by the Peres-Horodecki criterion, i.e. only 
“free” entanglement can be created in the PD process. 

Finally, we show that for a physically well-motivated 
family of /3e(0) states, generation of nonzero QEE at time 
t occurs iff the reduced state of E changes: pE{t) pe{0)■ 
Thus, measurement of change of any environmental ob- 
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servable can be a witness of creation of QEE. This is an 
interesting example how the common lore of the impos¬ 
sibility of bipartite entanglement detection by measuring 
one party only, turns out not to hold when some prior 
knowledge on the initial state and the character of inter¬ 
system coupling is available. 

We begin with the most general form of the Hamilto¬ 
nian of Q and E which describes the PD case: 

H = Hq + He + \0){0\(^Vo + \1){1\(^Vi . ( 1 ) 

The first term of the Hamiltonian describes the qubit and 
is given by Hq = X]i=o i second describes 

the environment, while the remaining terms describe the 
qubit-environment interaction with the qubit states writ¬ 
ten on the left side of each term (the environment oper¬ 
ators Vb and Vi are arbitrary). 

Since entanglement between two subsystems is unaf¬ 
fected by local unitary operations on either of the sub¬ 
systems [l3 - [l3| , the full qubit-environment evolution op¬ 
erator U{t) = exjp{—iHt) resulting from the Hamiltonian 
0 may be transformed into 

U{t) = = |0)(0|®lE+|l)(l|®w(i) , (2) 

where Hi = He + Vi and we have defined the operator 

w{t) = ejip{iHot) exp{—iHit) . (3) 

Note that while Hq commutes with all the other terms 
in H, this is not necessarily the case with He- 

To find out what the conditions for the generation of 
QEE via pure dephasing are, let us study the joint state 
of a qubit and an environment which are initially in a 
product state (t(0) = /5q(0) 0 /5e(0) and evolve accord¬ 
ing to the operator ([2]). The qubit is initially in a pure 
state Itjj) = a|0) -I- 6|1), with |ap -I- |5p = 1 and a, 6 ^ 0 
(a superposition is needed for dephasing to occur), so 
the density matrix pq(0) = |t/’)(V'| (the issue of QEE 
becomes more complicated when initial purity of Q is 
not maximal [^, and we exclude this case here). At 
this stage we impose no restrictions on the initial den¬ 
sity matrix of the environment and write it in terms of 
its eigenstates, pe(0) = The time-evolved 

qubit-environment density matrix in the evolution pic¬ 
ture described by the operator ([2]) takes the form 


|aPE„Cn|n)(n| ab*Y.n<^n\n){n'{t)\ \ 

- \ a*bY.^Cr.\n'{t)){n\ |5p c„|n'(t))(n'(t)| ) ’ 

( 4 ) 

where the matrix is written in the basis of the eigenstates 
of Q, and \n'{t)) = w{t)\n) with w{t) given by Eq. ([3]). 

Eirstly, let us study the simplest situation when the 
initial state of the environment is completely mixed 
(i.e. pe(0) = 1/fV where N is the dimension of the en¬ 
vironment). Using Eqs. ([2]) and (j4]) we arrive at 



|apl ab*lw^t) \ 
a*bw{t)l \b\^l ) • 


( 5 ) 


Now we use the complete mixedness of E: the unity can 
be written as a sum of projectors over any basis set, in¬ 
cluding the basis spanned by eigenvectors of a Hermi- 
tian operator h{t) defined by w{t) = exp[—ih{t)]. With 
h{t) \ k{t)) = Xk{t) \ k{t)) (note that the eigenvectors and 
eigenvalues evolve with time), we have 


a{t) 



|a|2 ab*e^^Vt) 
a*be-^^Vt) 


(g) \k{t)){k{t)\ 


which is by definition separable 0 and no QEE is cre¬ 
ated, even though the whole system can evolve in a com¬ 
plicated manner. 

In the general case, when /5e(0) is arbitrary, the 
question of QEE cannot be resolved so simply. The 
Peres-Horodecki criterion (PHC) [mu of negativity of 
partially-transposed a{t) detects some, but not all entan¬ 
gled states - the bound-entangled states Bin , which 
have a positive partial transpose, are much harder to de¬ 
tect, even for a constrained case of a system consisting 
of a qubit and an iV-dimensional environment 0 . We 
proceed now to check for the existence of negative par¬ 
tial transpose QEE. Later it will turn out that states 
not shown to be entangled by PHC are in fact separable, 
i.e. the bound QEE is never created in the process of pure 
dephasing of the qubit. 

The positive semidefinite matrix has only nonnegative 
eigenvalues, or equivalently all of its principal minors are 
non-negative. Checking for the second condition turns 
out to be manageable for the partial transpose of the 
density matrix from Eq. 0. Eirst we treat the case of 
a{t) which is full-rank (i.e. all c„ ^ 0). A class of prin¬ 
cipal minors can be obtained by symmetrically crossing 
out {N — 1) rows and columns in such a way that only 
one row and column containing a diagonal density matrix 
element proportional to |5p is left. They are given by 


/ |apco 


Mi = det 


a*bciy*Q \ 


0 ••• |apc„ a*bciy*^_-y 

\ab*Ciy^o ■■■ ab*c^yiN-l \b\'^J2n^n\yni\'^ J 


( 7 ) 

where i = 0,1,..., TV — 1 and yni = (n|w^(f)|i). A simple 
calculation leads to 


Mi = |ap^|6p X! ■ 

\ k / n=0 ^ 71 / 

( 8 ) 

When all Ci are equal [pe(0) oc 1 and entanglement 
is not generated as shown above], all Mi — 0 (this is 
obtained after noticing that both \yni\'^ = 1 and 
l?/mP = !)■ Let us assume now that not all c„ are 
the same. We choose then the minor Mi correspond¬ 
ing to the largest q, so that for all j we have Ci > Cj 
and at least for one j the inequality is sharp. We 
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have then Li = Yjn^ri\yni? < while 

2 

Pi = > CiX]„ly™P = Ci, in which we as¬ 

sumed that \yij\ ^ 1 for j corresponding to at least one 
Cj for which Cj < Ci. In such a case we see that Li < Pi 
and consequently Mi < 0 signifying the presence of non¬ 
bound (free) entanglement between Q and E. 

In the opposite situation, when yij = 0 for each j corre¬ 
sponding to Cj < Ci we repeat the same reasoning for the 
Mfc minor corresponding to the second-largest Ck- Then 
either we discover that the state is entangled, or y^j = 0 
for all j such that Cj < Ck- If iterating this procedure fails 
to sense entanglement (i.e. if no minors are shown to be 
necessarily negative), then it means that for any pair of 
i and j with Ci > Cj we have yij = 0, which is equiva¬ 
lent to (*| |j) = (j| w \i) = 0, so the evolution due to 

qubit-environment coupling does not lead out of the sub¬ 
space corresponding to Ci. In fact, due to unitarity of w, 
the non-negativity of all the Mi minors is equivalent to 
w having no matrix elements between states correspond¬ 
ing to different occupations c„. The evolution for which 
\n) —>■ |n) is a special case of this type when the sub¬ 

space in which the evolution takes place is limited to the 
single state \n). 

Let us show now that any evolution of this type is 
nonentangling. If we denote the subspaces by an addi¬ 
tional index s in such a way that any state jug) belongs 
to that subspace and all factors = Cg are equal, the 
qubit-environment density matrix from Eq. 0 can be 
rewritten as 


d{t) = ^Cg 

S 




(9) 

Since the subspaces do not overlap, each jug) (ug | is a 
unity in its own subspace (hence, also l^s(^))(^s(0l 
is the same unity), so the sum over s is a sum of den¬ 
sity matrices (up to a normalization) in which the re¬ 
stricted environment has an initial density matrix pro¬ 
portional to unity. It has been shown previously that 
such density matrices remain separable during PD evo¬ 
lution. Hence each of the matrices can be written in the 
form Qp^Ppsii) ® Rp,{t) and we get 


o-(t) = '^Cs'^qp^Pp,{t)^Rps(.t) = y^,Pkh{t)^Rk{t), 

s ps k 

( 10 ) 

with pk = CsQp^ and k numbering all the terms that are 
part of the summation (meaning that a single index k 
corresponds to a unique combination of the indices s and 
Ps). This shows that given the above conditions, the den¬ 
sity matrix a{t) is separable. The result concurs with the 
fact that entanglement between an N dimensional and an 
M dimensional system will always be free entanglement, 
if the number of eigenvalues of the joint density matrix 
is smaller or equal to max(fV, M) [20|. 

To complete the analysis of QEE generation via PD 


processes, let us study the situation when pe( 0) has 
eigenvalues equal to zero (it is not full rank). Obviously, 
if all states \p) corresponding to Cp = 0 are decoupled 
from all the other states, meaning that for all n for which 
Cn^O and for all p for which Cp = 0, (n|u)'f(f)|p) = 0, the 
Hilbert space of E can be reduced in such a way that it no 
longer contains the set of states {|p)}. The above proce¬ 
dure for finding entanglement can then be used for such a 
reduced environment and the same conclusions for when 
the evolution generates QEE hold. Let us now study the 
situation when there exists one state \p) for which Cp = 0 
but at least one {n\'w'^{t)\p) = ynp ^ 0 exists. The minors 
from Eq. ([5]) are now given by 


M, = -\a\^^\b\^ 



\yip\ 1 


( 11 ) 


for i p and Mp = 0. Hence, Mi corresponding to the 
non-zero value of yip is negative and free entanglement is 
generated. 

If there exist K > 2 environment states \p) for which 
Cp = 0, with at least one ynp ^ 0 for each, all the Mi given 
by Eq. ([5]) are equal to zero and they are no longer a good 
subset of minors for the study of QEE generation, yet 
there exists K new subsets, which can be used to probe 
QEE (the subsets are equivalent in the sense that each 
is sufficient to detect entanglement). These are obtained 
by symmetrically crossing out the same K — 1 rows and 
columns from each of the matrices given by Eq. 0 in 
such a way that only one diagonal element equal to zero 
•is left. If the state corresponding to this diagonal element 
is labeled as |r), the minors are given by 





c}\yv. 


( 12 ) 


As in the case of a single Cp = 0, the minor Mi cor¬ 
responding to the non-zero value of yir is negative and 
QEE is generated. 

This closes the part of the article where we character¬ 
ize the generation of QEE as a result of the general pure 
dephasing Hamiltonian O- We have shown that the ap¬ 
pearance of QEE at a given time depends only on the 
initial state of the environment, /5e(0), and on the form 
of the operator w{t) given by Eq. ([3]). Furthermore the 
qubit-environment state is separable at time t iff for all 
eigenstates |n) and |j) of pe{0) corresponding to eigen¬ 
values Ci^Cm we have {n\w{t) \i) = 0. This statement is 
equivalent to the statement that the qubit-environment 
state is separable at time t iff 

[pE(0),u;(t)] = 0 . (13) 

This is the main technical result of the first part of the 
paper. 
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The above criterion of course confirms the lack of 
QEE when /5 e(0)oc 1, since the unit operator commutes 
with any w{t). In the RU case the assumptions that 
[He , Vi] = 0 and pe = /{He) immediately lead to Eq. (IT^ , 
with w(t)ln} = (this is the case of a quasi-static 

E leading to inhomogeneous broadening of the energy 
splitting of Q). However, it is crucial to note that these 
examples do not exhaust the {pe(0), 'w(t)} pairs that sat¬ 
isfy Eq. (fT^ . This condition is fulfilled when pe( 0) has 
subspaces of equal c„, and w(t) operator is block-diagonal 
with respect to these subspaces. This means that qubit- 
induced evolution can be nontrivial within these sub¬ 
spaces in the absence of QEE. Eor example, one could 
take a system with pe(0) oc 1, and then by appropriate 
measurements on E post-select its states in such a way, 
that a subspace with altered Cn is singled out. If the 
qubit-induced evolution due to u)(t) does not couple this 
subspace with the rest of the E states, QEE remains ab¬ 
sent, while pe( 0) is nontrivial and nothing is assumed 
about commutation of He and Vi (except for the con¬ 
dition that there is a subspace closed with respect to 
evolution generated by w{t)). 

It is easy to see now that the transformed reduced den¬ 
sity matrix of E, 

PE(t) = TrQa-(t) = |ap/5E(0)-b|&pw(t)pE(0)-i5l'(t) , (14) 

remains constant for any a and &, iff the condition of 
separability from Eq. m holds: the lack of QEE is 
equivalent to pE(t) = Pe(0). After transforming back to 
the “laboratory frame” using Eq. m we see that when 
there is no QEE, the actual reduced state of E, pE{t) = 
exp(—i/foi)PE(0) exp(fifoO) remain unchanged iff 

[iJo,PE(0)] = 0. Our second main result is thus the fol¬ 
lowing: if this commutation condition holds, observation 
of any change in the reduced state of E during the evo¬ 
lution of the coupled qubit-environment system signifies 
the presence of QEE (note that all the above derivations 
could be repeated with Hi replacing Hq in Eq. m and 
the two Hi exchanging places in Eq. leading to anal¬ 
ogous condition involving Hi). In this case generation 
of entanglement can be verified by measurements per¬ 
formed only on one of the entangled subsystems (namely 
the environment), which is much more convenient than 
any viable joint qubit-environment measurements. Let 
us stress that in general (even when entanglement be¬ 
tween two qubits is studied), without any prior informa¬ 
tion about the initial state and the Hamiltonian, prob¬ 
ing the existence of entanglement via measurements on 
a single subsystem is not possible. Furthermore, in the 
discussed case the detection of any change in the state 
of E is equivalent to the detection of QEE. 

The constraint imposed on /5e(0) that enables the 
study of QEE via measurements on E only is not a very 
restrictive one. It means that Hi (with j = 0 or 1) can 
only have off-diagonal elements [when it is written in the 
eigenbasis of pe(0)] between the states corresponding to 


equal occupations c„. This occurs, for example, when 
/5e(0) is a function of one of Hi. The subspaces of equal 
Cn are then subspaces of equal energy. The most natural 
physical situation when this occurs is when the environ¬ 
ment thermalizes according to Hi. This can happen when 
one of the qubit-environment interaction terms, say Vq, 
vanishes (so that Hq = He), and E is allowed to reach 
thermal equilibrium in the absence of the qubit. This is 
the case for quantum dot excitonic qubit (with |0) corre¬ 
sponding to the optical ground state) coupled to a bath 
of phonons ^ [2^. Interestingly, in qubits based on 
the nitrogen-vacancy center in diamond (23| it is possible 
to realize both Ho = 0 and ^ 0 cases (by appropriately 
choosing two qubit levels out of three available states of 
electronic complex with total spin 1), and observe dis¬ 
tinct Q decoherence behavior in the two of them 24, 2^ . 
Alternatively, the qubit can be initialized in its lowest- 
energy state for a time long enough for E thermalize in 
its presence [i^, then the state of the environment is a 
function of the Hamiltonian Hq without the additional 
constraint of Hq = 0 (be. it reaches the thermal state 
defined by Hq before the qubit is rotated into a superpo¬ 
sition state and its dephasing begins). One should note 
that when thermal states of E are considered, the amount 
of generated QEE (presumably correlated with the mag¬ 
nitude of change of chosen environmental observable) is 
expected to be proportional to j3, the inverse tempera¬ 
ture of E (since for /3 = 0 there is no QEE). This could be 
useful for thermometry of an environment for which A:bT 
is large compared to typical energy scale of Hi, e.g. for a 
nuclear bath coupled to a spin qubit 27. 2^. 

In conclusion, we have studied the generation of qubit- 
environment entanglement (QEE) via pure dephasing 
processes and specified the only three classes of situations 
(specified by the initial state of E and relevant evolution 
operator which is derived from the Hamiltonian) when 
QEE will not be generated. These are, the case when the 
initial density matrix of the environment is proportional 
to unity, the case when the relevant evolution operator 
cannot change the occupation of any of the eigenstates of 
the density matrix (the random unitary evolution due to 
a quasi-static bath, essentially inhomogenoeous broad¬ 
ening of Q splitting in the context of pure dephasing), 
and a non-trivial mixture of the two cases which allows 
dynamical evolution within closed subspaces of equal oc¬ 
cupation. Furthermore, we have shown that restricting 
the initial E states to a class which is very common in 
any realistic qubit-environment setup, enables the use of 
a very powerful tool to measure QEE, since the state of 
E will remain static throughout the evolution if no QEE 
is generated. Hence, the detection of any change of the 
state of E is then equivalent to the detection of entangle¬ 
ment. 

This work is supported by funds of Polish Na¬ 
tional Science Centre (NCN) under decision no. DEC- 
2012/07/B/ST3/03616. L. C. thanks Igor Bragar for as- 
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sistance with numerical testing of hypotheses at an early 
stage of reseach, and Cezary Sliwa for useful comments. 


APPENDIX 

An example of evolution leading to entanglement of 
a qubit with maximally mixed environment 


Let us give a very simple example of a unitary qubit- 
environment operation (generated by a Hamiltonian 
more general than the one correspoding to pure dephas¬ 
ing of the qubit), which leads to creation of QEE at some 
time t after initialization of a pure state of Q and a com¬ 
pletely mixed state of E. We take a two-dimensional envi¬ 
ronment and take a unitary qubit-environment operation 
Uc that transforms a given qubit state together with 
the states |0), |1) of a two-dimensional environment as 
follows: 


ucm = joi), 

(15) 

UcIV'l) = (|00) + |ll))/y2 , 

(16) 

Uc\^±0) = im - \n))/V2, 

(17) 

Uc\^±l) = |10) , 

(18) 


where \ij) denotes the state \i)Q 0 |j)_E, and \tp±) is or¬ 
thogonal to \tj;). In the pure dephasing (PD) case the 
evolution is such that there exists a unitary operation on 
the qubit, that brings Uc to block-diagonal form with 
respect to states of Q. It is straightforward to check that 
such an unitary operation does not exist in the above- 
described case. Under the action of Uc an initial state of 
the form |^/')('!/'| 11e/2 is transformed into the state (in 

the qubit-environment basis {|00), |01), |10), |ll)}) 


Uc|V’)(V'l ® = 


/ 1/4 0 0 1/4 \ 

0 1/20 0 
0 0 0 0 
\ 1/4 0 0 1/4 y 


(19) 


which is entangled (the concurrence of this state is equal 
to 1/2). We have thus shown that there exist such initial 
states \^p) of the qubit, and such Hamiltonians of the 
whole system, that lead to creation of entanglement at 
some time t after initialization of the qubit and the two- 
dimensional environment in a separable state, with the 
environment being initially maximally mixed. 
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